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Abstract. A stress exists in solids surfaces, similarly to liquids, also if the 
underlying bulk material is stress-free. This paper investigates the surface stress 
effect on the measured value of the Si lattice parameter used to determine the 
Avogadro constant by counting Si atoms. An elastic- film model has been used 
to provide a surface load in a finite element analysis of the lattice strain of the 
x-ray interferometer crystal used to measure the lattice parameter. Eventually, 
an experiment is proposed to work a lattice parameter measurement out so that 
there is a visible effect of the surface stress. 
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1. Introduction 

The Avogadro constant, A^a, is determined by dividing the molar volume of a silicon 
single-crystal highly enriched with ^^Si by its unit cell volume, obtained from the 
measured value of the lattice parameter [H [21 E] • The lattice parameter is measured 
by means of combined x-ray and optical interferometry. Recently, the effect of surface 
stress was brought to our attention as a possible error source in our measurement of 
the lattice parameter [3]. Indeed, this issue was already taken into account but never 
definitely investigated [5]. An intrinsic surface-stress exists in solid, in a similar way 
as it occurs in liquids, which was first modelled by Gibbs [6] and Shuttleworth [7], 
whose works were seminal in surface thermodynamics. Though it is widely reported 
in review papers and textbooks on surface science [HI [9] , their theory is still a matter 
of debate [TU] . A continuous mechanics theory of solid surfaces has been developed by 
Gurtin and Murdoch [HI [12], which is summarized in p!^. When the surface relaxes, 
it compresses or pulls the underlying crystal; since the resulting lattice strain depends 
on the crystal size and geometry, the lattice parameter of an x-ray interferometer is 
different from that of the Si spheres used to determine the molar volume. 

After summarizing the operation of a combined x-ray and optical interferometer, 
this paper investigates the effect of surface stress on the interferometer operation. An 
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Figure 1. Sketch of a combined x-ray and optical interferometer. The analyser 
crystal moves in a direction parallel to the laser beam, its displacement is measured 
by the optical interferometer, and each passing diffracting plane is counted. Both 
the x-ray and optical interference patterns are imaged into position-sensitive 
detectors to sense lattice strain and parasitic rotations of the analyser. 



elastic membrane model has been used to provide a surface load in a finite element 
analysis of the lattice strain. Particular emphasis was given to the operation of the 
^^Si interferometer used to determine A^a- Eventually, a measurement of the lattice 
parameter with a variable-thickness interferometer is proposed in order to evidence 
the effect of surface stress on the underlying crystal. 

2. Combined x-ray and optical interferometry 

A review of the lattice parameter measurement can be found in [14j. As shown in 
Fig. [U an x-ray interferometer consists of three Si crystal slabs so cut that the {220} 
planes are orthogonal to the crystal surfaces. X-rays from a 17 keV Mo Ka source 
are split by the first crystal and recombined, via two transmission crystals, by the 
third, called the analyser. When the analyser is moved in a direction orthogonal to 
the {220} planes, a periodic intensity-variation of the transmitted and diffracted x- 
rays is observed, the period being the diffracting-plane spacing, about 192 pm. The 
analyser embeds front and rear mirrors, so that its displacement and rotations can be 
measured by an optical interferometer. The measurement equation is ^220 = ^A/ (2n), 
where ^220 ^ 192 pm is the spacing of the {220} planes and n is the number of x-ray 
fringes observed in a displacement of m optical fringes having period A/2 ^ 266 nm. 
To ensure the interferometer calibration, the laser source operates in single-mode 
and its frequency is stabilized against that of a transition of the ^^^l2 molecule. To 
eliminate the adverse influence of the refractive index of air and to achieve millikelvin 
temperature uniformity and stability, the experiment is carried out in a thermovacuum 
chamber. Continuous developments led the measurement sensitivity and accuracy to 
approach 10~^(i22o; the most recent results are given in [3]. 

3. Surface stress 

Similarly to liquid surfaces, solid surfaces exhibit a stress also if the bulk material is 
stress-free. This is the macroscopic consequence of atomic affects: roughly speaking, 
the atoms on the surface lose half of the bonds they would have in the bulk; therefore, 
crystal surfaces are disturbances in the bonding of atoms. 



3 



Table 1. Calculated values of the principal components, (7|| and cr^, and mean 
value, (7001, of the surface stress of the symmetrically dimerized Si(OOl) (2 x 1) 
surface. The || and _L directions are parallel and orthogonal to the dimer rows; 
positive and negative values indicate tensile and compressive stresses. 
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Surface stress was first generalized to solids by Gibbs and Shuttleworth |6| i7j as 
the amount of reversible work per unit area needed to stretch a pre-existing surface. 
Since a solid surface may change not only by exchange of matter with the bulk, but 
also by elastic deformation, the concept of surface stress must be derived from the 
molecular dynamics of the surface [151 US] • 

Continuum mechanics provides the framework to define the surface properties 
and the relationships among them. The elasticity theory of solid surfaces to which we 
refer was developed by Gurtin and Murdoch [HI [12] . A summary of the linearised 
theory - which is valid when both the surface and bulk strains are small - is given 
in [13]. The key concept is that, for a surface attached to the underlying solid, a 
stress- free reference does not exist; since the bonds between the surface atoms and 
between the bulk atoms are different, the surface of a stress- free crystal has residual 
strain, eg^, and stress, ctq^. 

If the surface tends to shrink (expand) with respect to the bulk, the surface stress 
is positive (negative) and it is said to be tensile (compressive). Since a fully relaxed 
surface has no normal stress, surface strain and stress are second rank tensors, which 
are characterized by two principal orthogonal-components. If the surface has more 
than a 2-fold symmetry, the surface stress and strain tensors are isotropic. In this 
case, = eoSai3 and the free surface-energy per unit area is [13j 

7 = 70 + 2(/i + X)eoeaa + l^^a(3^a/3 + XCaa^/S /3 f'^, (1) 

where 6^(3 are the strain component in the plane tangent to the surface, 70 is the 
surface energy for a stress-free bulk, fi and A are elastic constants, and the summation 
over the repeated indices - a, /3 = 1, 2 - is implied. The definition of surface stress is 
o'a/s = dj/dcap] hence, by taking the derivatives, 

0'a(3 = CroSa^ + 2/16^/3 + ^i^aa + e/3/3)^^a/3/2, (2) 

where ctq = 2(/i + A)eo is the intrinsic surface stress. 

The determination of the surface stress in solids is a challenge; the measurement 
of the lattice parameter in small particles has been the main method Given the 
experimental difficulties and since the results are affected by the sample preparation 
and adsorbed impurities, the values obtained from ab-initio and molecular dynamics 
calculations seem more reliable. The results for the symmetric (2x1) reconstruction of 
the Si(OOl) surface are summarized in table [TJ the reconstructed surface is formed by 
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domains of anisotropic stress, but, for symmetry reasons, the mean stress is isotropic. 
The scatter of data suggests a reassuring smah surface stress. 

Because of the anisotropic nature of the crystal, surface stress depends on the 
surface orientation. A measurement of the polar dependence of the relative surface- 
stress at high temperature (1373 K) is given in [23j, which allows the surface stress 
of any surface to be determined from that of the (001) reference surface. In the case 
of our x-ray interferometer, the (110) surfaces are those of main interest; because of 
their 2-fold symmetry, two values are necessary to represent the surface stress of the 
two perpendicular [001] and [110] components. As reported in |23j, at 1373 K, the 
sought principal components are ctj^^ ~ 0.5<jooi and ctj^^ ~ 0.3crooi- Indicative aooi 
values are given in table [TJ 

The crystal's surfaces are normally covered by a few nanometre thick oxide- 
layer, which can further stress the surface. The stress changes of the Si(OOl) (2 x 1) 
surface during plasma oxidation were measured by using a thin silicon cantilever; the 
bending caused by the oxidation of one face was optically measured and related to the 
surfaces stress difference by elasticity theory [23]. Initially, the oxidation induces 
a compressive stress, followed by a tensile stress; eventually, the stress variation 
turns to compressive for oxide thickness greater than about 1.5 nm. The initial 
stress changes were explained in terms of the electron diffusion into the sample and 
the shrinkage caused by the oxygen-bridged dimer structure. The final compressive 
variation, A<Jooi ^ —0.1 N/m, was explained by the oxidation of Si atoms at deep sites; 
since the oxygen expands the Si-Si bond, oxidation is expected to cause a compressive 
stress. 

The surface relaxes the intrinsic stress by compressing or pulling the underlying 
crystal. For instance, the (isotropic) lattice strain of a sphere of radius R is 

2c^o _ 2ao 

^~ 4(/i + A) + 3Ki? 3KR' 

where K is the bulk modulus and the Young-Laplace approximation holds when the 
elastic constants can be neglected when compared to the KR product. By choosing 
(jQ = 1 N/m, we can use this formula - where 2R = 93 mm and K = 98 GPa 
- to estimate the lattice strain of the spheres used to determine A^a- The result 
is e = —0.07 nm/m, which is irrelevant for the A'a determination. However, if 
2R =1.2 mm - which is a value equal to the thickness of the crystal used to measure 
the lattice parameter - we obtain e = —6 nm/m, which is significant with respect to 
the uncertainty of the lattice-parameter measurement. 

4. Finite element analysis 

Figure [5] (left) shows the analyser crystal of the ^^Si x-ray interferometer used for the 
A'a measurement and its crystallographic orientation. The end mirrors are parallel to 
the diffracting planes; to avoid stress propagation, they are separated by vertical cuts. 
The lamella thickness, 1.2 mm, was chosen by a trade off between the increased lattice 
strain (due to the surface stress) of a thin lamella and the increased x-ray absorption 
of a thick lamella. 

To estimate the lattice parameter change caused by the surface stress, we used a 
commercial finite element analysis software [23]. As shown in Fig. [21 the (220) planes 
are orthogonal to the lamella faces which are parallel to (110) planes. Owing the large 
uncertainty of the residual stress values, though the (110) surfaces are anisotropic. 
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Figure 2. Left: sketch of the analyser crystal of the ^^Si interferometer. M: end 
mirrors, A: analyser lamella, B: base. The (5 x 5 x 0.1) mm^ reliefs that support 
the crystal are coloured in black. The lattice parameter was measured along the 
indicated line, a.t z = 14 mm. Right: surface-stress modelling by shear forces. 

the analyser surfaces were modelled as stretched by isotropic membranes attached to 
the underlying crystal lattice. Since we are only concerned with infinitesimal strains, 
though the surface stress aaf3 depends on the strain according to ([2j), we used an 
isotropic ctq = 1 N/m value. This surface stress value allows the analysis results to be 
easily scaled if a different value has to be considered and, hopefully, it is larger than the 
actual stress value. Hence, according the Young-Laplace equation, the surface applies 
the inward pressure P = 2/^(7o, where is the mean curvature. As schematically shown 
in Fig. [5] (right), when considering flat surfaces, the surface stress was modelled by 
shear forces per unit length applied orthogonally to the surface edges. 

As strains are inflnitesimal, we uses a linear model; the mesh, of about 10^ 
tetrahedral anisotropic elements, was the result of successive reflnements. The not 
null elements of the stiffness matrix are given table [5] in both the [100], [010], [001] 
and [110], [110], [001] crystal axes. As regards the model parameters, the Si density 
and gravitational acceleration are 2330 kg/m^ and 9.81 m/s^, respectively. 

The analyser rests on a silicon plate; it stands on three (5 x 5 x 0.1) mm^ reliefs 
carved out of the crystal base by chemical etching. In order to avoid stresses, it is 
held in position by a thin layer of high- viscosity silicone oil. Since the oil applies only 
viscous forces, the contact with the support has been simulated by setting to zero the 
displacement of the contact nodes. 



Table 2. Not null elements of the stiffness matrix for a silicon crystal. The 



measurement unit is GPa. 


[100], [010], [001] axes 

Cll,C22,C33 C12,C13,C23 C44, C55 


C66 


165.7 63.9 79.6 


[110], [110], [001] axes 

'^ll''^22 ^12 ^13' ^23 ^^33 ^44' ^^55 


^66 


194.4 35.2 63.9 165.7 79.6 


50.9 



Figure 3. Left: exx component of the diffracting-plane strain for the optimal 
location of the supports, where blue is —10 nm/m and red is 10 nm/m. Right: 
qualitative magnified image of the diffracting planes. A realignment of the 
analyser has been simulated by subtracting the average displacement and tilt. 
The displacement magnification factors is shown by the 1 nm bar. 



4.1. Results 

We calculated the strain Cxx = dxUx-, where Ux is the horizontal component of the 
displacement vector. Both the effects of self- weight and surface strain have been 
considered; the supports were so chosen as to minimize bending or sagging The 
three dimensional map of the strain distribution in the analyser and a magnified image 
of the strained diffracting planes are shown in Fig. [3l 

The Cxx strain is trivially related to the lattice parameter by A(i22o/<^220 = ^xx- 
Hence, the relative variations of the 6^220 values along measurement lines at different 
heights are identical to the Cxx values calculated along the same lines; they are shown in 
Fig. m (left). The figure indicates that, in the measurement horizontal strip a couple of 
millimetre wide at a z = 14 mm, the lattice is uniformly strained by about —6 nm/m, 
as estimated by the previous "spherical-cow" model. In addition to the horizontal 
strain, there is a vertical strain gradient of about 3 nm m~^ cm~^, the bottom 
part being less strained, which could be the subject matter of future experimental 
verifications. 

Figured] (right) updates the corresponding figure 6 of [2B], which was used to 




Figure 4. Left: horizontal ci220 variations in y = plane, at different distances 
from the analyser bottom-surface. Right: Relative d220 variations along the 
z = 14 mm measurement baseline (figure[2]). Different support-to-mirror distances 
have been considered, as indicated in the graph. 
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Figure 5. Sketch of the variable-thickness analyser; vertical and horizontal cuts 
separate the end, thin, flag-like lamellae from the mirrors and base. The expected 
d220 variations are calculated along the indicated lines, at z = 9 mm. 

assess the sole effect of self weigh in the lattice parameter measurement. It shows 
that, if the analyser surfaces are loaded by a tensile stress of 1 N/m, the 6^220 value 
is about 6 x 10~^(i220 smaller than its value in an unstrained crystal. The figure also 
shows that the lattice parameter profile given in |3j cannot deliver clues about the 
surface contribution to the measured ^220 value because the only effect of the surface 
stress is to translate all the curves downwards by 6 x 10~^(i220- 

Since the uncertainty associated to the measured ^220 value is 3.2 x 10~^(i220, 
the lattice strain induced by a 1 N/m stress impairs the measurement accuracy. The 
results indicate also that we can cope effectively the surface stress by a numerical 
determination of its effect. However, owing the large uncertainty of the sign and 
value of the surface stress, we don't yet propose to reconsider the corrections and 
error-budget contributions given in [3] . Further investigations are under way to assess 
and to quantify the effect of surface stress, if any. In the next section we describe 
an analyser design to work a lattice parameter measurement out so that there is a 
measurable effect of the surface stress. 

5. Variable thickness interferometer 

A way to evidence if the surface stress strains the analyser crystal up to an extent 
which is significant to the lattice parameter measurement is to compare the results 
of measurements carried out in crystals having a large thickness difference. Since the 
induced ^220 changes are small, to avoid they are masked by lattice imperfections, 
measurements must be carried out in the same crystal. Consequently, as shown in 
Fig. [5l we designed an x-ray interferometer having a variable analyser thickness. The 
minimum 0.4 mm thickness is set by machining capabilities; the maximum 1.5 mm 
thickness is set by x-ray absorption. 

After the optimization of the finite element analysis and the location of the 
supports so as to minimize the self- weight bending, we obtained the strain distribution 
shown in Fig. [6] (left); a magnified image of the strained {220} planes is shown in Fig. [6] 
(right). The expected ^220 profile along the horizontal lines at z = 9 mm height is 
shown in Fig. [3 (left). As shown in Fig. [3 (right), the surface stress and self- weight 
have opposite effects. 
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Figure 6. Left: e^x component of the diffracting-plane strain for the variable- 
thickness analyser, where blue is —10 nm/m and red is 10 nm/m. Right: 
qualitative magnified image of the diffracting planes. A realignment of the 
analyser has been simulated by subtracting the average displacement and tilt. 
The displacement magnification factors is shown by the 1 nm bar. 

The surface stress will be estimated from the (^220 gap between the thick and thin 
analyser parts. If the surface is loaded by a tensile stress of 1 N/m, ^220 decreases by 
about 15 X 10~^(i220 from the central (thick) to the outer (thin) parts. This variation 
is large enough to be detected. With a high crystal perfection, the present detection 
limit of 3x 10~^(i220 corresponds to a minimum detectable stress of 0.2 N/m. To ensure 
that self-weight bending does not change in a significant way the ^220 gap between 
the thick and thin analyser parts, the height of the analyser base was increased to 
19.8 mm. Figure [H] shows that, when the analyser supports are randomly displaced 
(m. X — y plane) from their optimal positions up to ±1 cm, the ^220 gap between the 
thick and thin analyser-parts changes less than what can be experimentally detected. 

6. Conclusions 

The influence of the intrinsic surface-stress on the operation of the x-ray interferometer 
utilized to measure the lattice parameter of the ^^Si crystal used to determine the 
Avogadro constant has been investigated by finite element analysis. The crystal 
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Figure 7. Horizontal ^220 variations in y = plane, along the 2 = 9 mm 
measurement lines (figure [Sj. The supports are optimally located so as to 
minimize the self- weight bending. Left: Both the self- weight and surface stress 
are considered. Right: Only the self- weight is considered. 
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Figure 8. Difference between the horizontal d220 variations in y = plane, along 
the z = 9 mm measurement lines, when the analyser supports are displaced (in 
X — y plane) from their optimal positions up to ±1 cm. 

surfaces, modelled as membranes having a tensile stress of 1 N/m, relax by compressing 
the underlying crystal. Though the induced strain is greater than the uncertainty 
associated to the measured value of the lattice parameter, since the stress magnitude 
is highly uncertain and the chosen 1 N/m value could be pessimistic, we do not propose 
a correction of the published value. To verify if the surface stress contributes or does 
not contribute to the lattice parameter measurement, a variable-thickness design of 
the x-ray interferometer has been investigated to induce detectable strain variations 
in the same interferometer crystal. The realization of such an interferometer is under 
way in collaboration at the Physikalisch Technische Bundesanstalt in collaboration 
with the Leibniz-Institut fiir Oberflachenmodifizierung. 
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